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Measurements provide a novel mechanism for generating the entanglement resource necessary for 
performing scalable quantum computation. Recently, we proposed a method for performing parity 
measurements in a coupled quantum dot system 1 . In this paper we generalise this scheme and 
perform a comprehensive analytic and numerical study of environmental factors. We calculate the 
effects of possible error sources including non-ideal photon detectors, ineffective spin-selective excita- 
tion and dot distinguishability (both spatial and spectral). Furthermore, we present an experimental 
approach for verifying the success of the parity measurement. 

PACS numbers: 



I. INTRODUCTION 

Quantum dots are semiconductor hetero-structures 
which exhibit strong electron and hole confinement. This 
leads to a highly discrete level structure, giving rise to 
many interesting nanotechnological applications, such as 
quantum information processing (QIP). Such a discrete 
level structure enables us to identify well-defined effec- 
tive two- level systems (qubits), which we use to encode 
our quantum information. A natural qubit is the spin 
of an excess electron within a quantum dot, which typ- 
ically exhibits long lifetime (up to milliseconds^) and 
coherence times (up to microsecon dP). 

For universal quantum computation (QC) we must be 
able to perform any one qubit gate, plus an entangling 
two qubit gate. Typically we utilise the natural couplings 
between the spins to generate our two qubit gate. Exper- 
imentally these two qubit gates would be achieved as pe- 
riods of free evolution of the interacting system followed 
by periods of controlled single qubit rotations. For ex- 
ample, early proposals 5 included using the the exchange 
interaction H = J Si • 62 to provide the necessary two 
qubit interaction. However, precise manipulation of the 
interactions between the spins is needed for such a scheme 
to work, which might be difficult to achieve in practice. 

An alternative method, which we discuss here, is to 
utilise the power of measurements to generate inter- 
actions between the qubits. Such measurement-based 
ideas initially appeared in the context of linear optical 
approaches to quantum computation^. In 2000, Knill, 
Laflamme and Milburn proposed a method for generat- 
ing entanglement using only linear optical elements, an- 
cilla photons and measurements on single photons. A 
number of proposals have since been presented for gen- 
erating entanglement using optics in solid-state systems. 
For example, Refs. [71 and [8] propose using single pho- 
ton interference to entangle spatially separated matter 
qubits. 

If we could implement a measurement-based scheme 
with the electron spin we could do away with the 



exquisite control required for the conventional ap- 
proaches discussed above. However, there exists a no-go 
theorem which states that it is not possible to achieve 
an exponential speed-up over classical computation us- 
ing solely single electron Hamiltonians and single spin- 
measurements 9 . Recently, Beenakker et a/P^have shown 
that it is possible to lift this restriction if we look out- 
side the Hilbert space of a spin and exploit the charge 
degree of freedom. Charge and spin commute, and so we 
are able to make measurements on the charge without 
destroying any information that is contained in the spin 
degrees of freedom. Beenakker et al. proceed to show 
that partial-Bell state measurements (also known as par- 
ity measurements) on the spin states are sufficient to im- 
plement a CNOT gate, thus enabling universal QC. A 
number of subsequent papers have proposed specific im- 
plementations for these spin-parity measurements. These 
include a charge tunneling detection method^), and a 
charge fluctuation methocM 

Recently, we proposed a method for optically perform- 
ing parity measurements on a pair of coupled quantum 
dots. In this paper we generalise this scheme and perform 
a comprehensive study of environmental factors. We will 
begin, in section [IT] by outlining the system and the inter- 



actions present. In section III we will detail the various 
steps in the spin-parity measurement. We will present 
results for the operation of the measurement in section 
[TV] , and then in section |V| we will analyze the effect of 
various error mechanisms, including valence band mix- 
ing, and spatial and spectral distinguishability of the two 
dots. We then proceed to describe a method for verifying 
the success of our gate in section [VT] . Finally, we will give 
some concluding remarks. 
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FIG. 1: Energy level diagram for the coupled quantum dot 
structure. Resonant transitions are denoted by the solid lines, 
while non-resonant transitions by dotted lines. 



II. THE COUPLED QUANTUM DOT SYSTEM 

In order to develop our coupled model, we will first 
introduce the structure of a single quantum dot interact- 
ing with a laser field. We assume that the dot is n-doped 
such that the valence levels (VL) are completely filled 
and the only occupied conduction level (CL) is the low- 
est lying level. In what follows, we will only consider 
the top most filled valence states \J Z = ±3/2). The 
CL electron has a spin degree of freedom on which we 
encode our quantum information: |0) is encoded in the 
m z = — 1/2 state and |1) in m z = 1/2. The system is ir- 
radiated by a classical a + circularly polarised laser field, 
resonant with the VL-CL energy gap. An electron-hole 
pair (exciton) state can be created, if a photon's angu- 
lar momentum of +% can be absorbed; Pauli's exclusion 
principle means this is only possible when the qubit is 
in state \m z = -1/2) = |1). This effect is Tauli Block- 
ing', and it enables us to generate excitons conditioned 
on the state of the qubit electron. The combined qubit 
electron/exciton state is a trion, denoted by \X). 

When two such doped dots are placed close to each 
other, there are direct electron spin- spin couplings. How- 
ever, these are very weak: current experiments place their 
strengths at less than 1 /ieV^P. We take these interac- 
tions to be negligible in comparison to the excitonic inter- 
action that we will exploit. The two primary excitonic 
couplings that we consider are a static and a dynamic 
dipole-dipole coupling. The static interaction results in 
an energy shift Vxx of the double trion state and is 
due to the inter-dot exciton-exciton dipole interaction. 
The dynamic coupling, or resonant Foerster interaction 
Vp, results in an exciton transfer from one dot to the 
other that is mediation by a virtual photon. This pro- 
cess has been shown, to first-order, to be non-magnetic 
and thus conserves electron spirl^. Therefore, in this sys- 
tem the Foerster interaction only couples the \X1) and 
| IX) states. The Hamiltonian for the coupled QDs is: 



H = uj \X)(X\®I + ujoI®\X}(X\ + 

Vxx\XX){XX\ + V F (\lX)(Xl\+H.c.) 
+Slcosu;it(\l)(X\®i + i®\l)(X\ +H.c 



,(1) 



where H.c. denotes Hermit ian conjugate, uoq is the exci- 
ton creation energy for both dots (the dots are assumed 
to be identical), £1 is the time-independent laser coupling 
(assumed to be the same for both dots), and uo\ is the 
laser frequency. The energy difference between the |0) 
and |1) states is negligible on the exciton energy scale. 

The Hamiltonian may be decoupled into four sub- 
spaces with no interactions between them: Hqq = {|00)}, 
Hoi = {|01),|0X)} , H 10 = {|10),|X0)}, H n = 
{|11), \X1), \X1), \XX)}. Let us look more closely at the 
Hamiltonian for the last of these subspaces. We rewrite 
this in basis of the eigenstates when Q — 0, which are 

|ii), |v+> = ^(|ix) + |xi», |v._> = ^(|ix) - |xi» 

and | XX). The degeneracy of the |^_) and |^ + ) levels is 
lifted by the Foerster interaction, resulting in two states 
each containing a delocalized exciton. In this basis the 
Hamiltonian is, after re-introducing a finite Q: 



H u = (o; + V»|^ + )(^+| + (o; - V»|^_)(^_| 
+{2uv + Vxx)\XX){XX\ 

+ti cosc^(|ll)(^+| + \rl>+)(XX\ +ff.c.).(2) 

The only dipole allowed transitions in this subspace are 
between |11) and 1*0+), and between |^ + ) and \XX), with 
a coupling strength of ft = ^Vt. The level structure for 
all four subspaces is shown in Fig. [I] 



III. SPIN PARITY MEASUREMENT 

The parity measurement protocol consists of two steps: 
excitation followed by monitored relaxation. In the ex- 
citation we aim to transfer the population of the states 
from the computational basis into the excitonic levels in 
the odd-parity (|10) — 1 01) ) subspace, while retaining the 
population of the even-parity ( 1 00) — |11)) states in the 
ground levels. We can achieve this by exciting the cou- 
pled dots with a pulsed laser tuned to energy c^o- By 
referring to Fig. [T] and Eqs. [2] and [I] we can see that 
such a laser excites transitions within the 1 01) and 1 10) 
subspaces. The 1 1 1 ) state is not excited to first order if 



\v F \,\v xx \^\n |/2, 



(3) 



since the laser is off resonance with the excited levels in 
this space, and 1 00) is of course optically inactive. 
Suppose we have an initial state 

|^o) = «oo|00) + aoi|01) + aio|10) + an|ll). (4) 

After an excitation tt pulse, our state is = aoo|00) + 
aoi|OX) + aio|XO) + an|ll). We next enter a period 
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of monitoring the system for decay photons. Assuming 
perfect detection, the state of the system is projected 
into the odd parity subspace if a photon is detected, and 
into the even parity subspace if a photon is not detected. 
Importantly, we only distinguish between the two parity 
subspaces without distinguishing states within the same 
subspace. We can represent the action of the measure- 
ment in terms of projection operators for the two desired 
outcomes: 



that there is no information gained about the source of 
the photon that is emitted. There are many ways in 
which this condition can be compromised: for example 
spatial and spectral distinguishability of photon emis- 
sions from the different dots. Our system will also have 
imperfect detectors; in the next section we will begin to 
look at these potential sources of error. 



Po 
Pe 



\^odd){^odd\ 

eve 



|01)(01| + |10)(10|, 
= |00)(00| + |11)(11| 



(5) 



As mentioned earlier, we must not be able to distin- 
guish between states within the same subspace. There- 
fore it is important throughout the radiative relaxation 



IV. CONDITIONAL DYNAMICS WITH 
IMPERFECT DETECTORS 

To model the monitored radiati ve relax ation we use 
the quantum trajectories formalis m 115 * 16 * 17 !. The condi- 
tional master equation (CME) describing n monitored 
relaxation channels is: 



n 

dp c = -i[H,p c ]dt + ^2 [vj Tr (J[cj]pc)Pc + (1 - Vj)J[cj]pc ~ A[cj]p c }dt 



f J[cj]pc 
\Tr{J[ Cj ]p c ) 



:}dNj(t) (6) 



where p c is the density matrix of the system, H is the 
system Hamiltonian in the interaction picture, Cj is the 
Lindblad operator through which the system couples to 
the measurement channel j, J[cj) is the jump super- 
operator which projects out the component of the state 
that is consistent with a detection from channel j and 
is defined as J[cj\p c = Cjp c Cj. A[cj] is defined to be 
A[cj]p c = ^(cjCjPc + p c CjCj); r]j is the efficiency of the 
detector that monitors emission into channel j. dNj(t) is 
the classical stochastic increment taking the values {0, 1}, 
which denotes the number of photons detected in channel 
j in the interval t,t + dt. 

Between quantum jumps, when dN(t) = 0, Eq. [6] is 
equivalent to the linear, unnormalised, CME 



Parameters 


Value 


V F 


0.85 meV 


Vxx 


5 meV 




2 eV 


Vt 


0.1 meV 


rx 


1 ns 


T x 


4 fieV 



TABLE I: Table of relevant parameters for the coupled quan- 
tum dot system 



For an initial state Q, and using Eq. [7| we find that 
the probability pe that we are in the even parity subspace 
at a time t after excitation is: 



■p = - i[H , p) + {(1 - Vj)J[cj]P ~ Acj]P) (7) 

3 

where p c = p/Tr(p). 

Assuming that the condition ([3| is satisfied, there are 
no excitons in the even-subspace. Thus we need only 
consider one channel defined by the Lindblad operator 
c = vTx(|01)(0X| + |10)(X0|). Furthermore, the odd- 
space only contains a single excitation and thus we can 
model the dynamics in two steps: a period of continuous 
evolution followed by a potential single quantum jump 
due to the photon detection event. 

Suitabl e parameters are chosen and are displayed in 
the Table [iVl 8 :* 19 * 20 * 2 ^ 
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(8) 



In order to increase this probability, it makes sense to 
wait for long enough that we can be sure that if a photon 
has not yet been emitted, it is not likely to be emitted 
in future. This amounts to waiting for a time t ^> 1/Yx- 
Then the fidelity of projection into the even-subspace is: 



F E 



a 



li 



(1 - rj)(a 2 01 + aj ) + K + ai ± ) 



(9) 



When a photon is detected (and ignoring typically neg- 
ligible detector dark counts) the fidelity of for projection 
into the odd-subspace is Fo = 1 . 

The success of our two step parity measurement is 
strongly dependent on detector efficiency, and can be- 
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come quite poor for typical values of rj. However, by re- 
peating the spin-parity measurement one or more times, 
it is possible to obtain improved fidelities. On each round 
of the repeated measurement, we gain greater confidence 
that we have successfully projected into the even sub- 
space, rather than missing every emitted photon. To 
analyze this we write the effect of the spin parity mea- 
surement when no photon is measured in the quantum 
operation formalism^. The action of a general quantum 
operation can be written as: 



(10) 




0.4 0.6 0.8 

Efficiency r\ 



where the are the Krauss projection operators which 
must satisfy the normalistion condition ^ k E^Ek = I- 
For our particular example, a single operation of the spin 
parity measurement will yield: 



c E pE' 



■ c E E E pE E 



Tr[c E pE ] 



c E E E pE E ] 



(11) 



where the un-normalised Krauss operators, defined as 
Eo,e = Po,e = \^o,e)(^o,e\' are tne Projectors onto 
the odd and even parity subspaces respectively, and the 
coefficients are chosen in a self-consistent manner: ce = 1 
and Co — (1 — rj). The denominator in this expression 
ensures that we satisfy the normalisation condition. 

When the measurement is repeated r times, the overall 
quantum operation is given by: 



p - £{£{p)) = 



c r E pE^ + c r E E E pE E 



Tr[c r E p& + c E E EP E E 



t- 



(12) 



The fidelity of correctly projecting into the even subspace 
is therefore given by: 



Fl 



Tr[c E E E pE E ] 



Tr[<f EopEl 



(ago ^ ~ 2 



c E E E pE E ] 
n) 



l 10J 



) + {a 2 m + a\ 1 ) 



(13) 



For all non-zero detector efficiencies, in the limit r — > 
oo the term (1 — r)) r will tend to zero. We can therefore 
expect a unit fidelity for every input state in this limit. 

Let us now look at the average fidelity for all input 
states as a function of r. Owing to the normalisation 
condition Oq Q + Qqi + a io + a ii = 1, we write the four 
coefflcents in terms of four-dimensional hyperspherical 
polar coordinates: 



«00 


= sin c 


h 


a 01 


= COS ( 


h 


aio 


= sine; 


h 


an 


= sine; 


h 



(14) 



FIG. 2: The averaged fidelities Fe as 
efficiency. 



and the area element is given by: 

dA - sin 2 (j>i sin <p 2 d(j)\ . 
The resulting integral is: 



function of detector 



(15) 



d<pi sin 



2tt 



d(j) 2 sin 02 / dfaFz. (16) 



This integration is performed numerically and the re- 
sulting averaged fidelity as a function of r\ can be seen 
in Fig. |2j We clearly see a convergence, as r — > oo, of 
the average fidelity to unity for all non-zero detection ef- 
ficiencies. Thus, by simply repeating the measurement, 
we are able to overcome the inherent problems of loss 
detectors. 



V. CORRECTIONS TO MODEL 
A. Valence band mixing 

The Pauli blocking mechanism, crucial to the success 
of the excitation step, is valid only in the case of no 
light-heavy hole mixing, which is only true in very lim- 
ited casesp^. Usually the hole eigenstates are composed 
of mixtures of light hole (\J Z — ±1/2)) and heavy hole 
(\J Z = ±3/2)) states^. The mixing is characterised by a 
factor e. This results in two families of eigenstates: one 
with predominantly light hole character and one with 
predominantly heavy hole character. The latter tend to 
be the topmost valence levels and are^t 



\h+) = Vl~e 2 \J z = ±3/2) 
\h-) = v / T T ^ = -3/2) 



■€^ = -1/2) 

■e\J z = +l/2). (17) 



On applying the <r + polarised laser field, it is now possible 
to generate excitons from both the |0) = \m z = —1/2) 
and |1) = \m z = ±1/2) electron states. The resulting 
exciton-laser field coupling Hamiltonian is: 



H a+ = ca8[u L t](\l)(X-\ + e\0)(X+\+H.c.) 



(18) 
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FIG. 3: The effect of hole-mixing on the level structure of the 
coupled quantum dot system. The two Foerster interaction 
terms are Vf — Mhh,hh and Vf = 2eMih,hh- Once again, far 
off resonant transitions are denoted by dashed lines and near 
resonant transitions by solid lines. 
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where the trion levels are 
modified mixing angle is e = e and khhh are charac- 
teristic lengths associated with the overlap of the electron 
and the Ih, hh Bloch functions (see Ref. Hl)l . 

The mixing also induces a Foerster interaction that 
couple other single exciton levels: 



Hp 



M hhM (\0X + )(X + 
2MiH hh e 



V3 



0X_)(X+1| + |1X+)(X_0|) + h.c. 

(19) 



where Mij are the matrix elements for the transitions 
induced by the Foerster interaction and i and j denote 
the different initial and final hole states respectively. 
The Hamiltonian for the coupled quantum dots is now: 

H = ^ o (|X + )(X + | + |X_)(X_|)0/ + 
<j I®(\X + )(X + \ + \X-)(X-\) 

+H F + H a+ + Vxx\^)(HW 
v,ve{X-,x + } 

With the inclusion of hole-mixing each decoupled sub- 
space will now have two single excitonic levels and a biex- 
citon level, and this is illustrated in Fig. [3j We now dis- 
cuss whether we can still perform the parity projection 
in this more complex situation. 

First, we must ensure that transitions to excitonic lev- 
els in the even subspace remain suppressed. Condition 
[3] is still valid to for suppressing transitions from |11). 
Meanwhile, the couplings to the excitonic levels from the 
1 00) state are reduced by the mixing factor e, and so con- 
dition [3] is also sufficient to suppress these transitions. 

Second, transitions within the odd subspace must only 
occur between the zero and single exciton levels. Our 
measurement replies on the detection of a single photon: 



on detection of the photon any population in the biexci- 
ton level will be projected into the single excitonic levels, 
and this must be avoided. 

The dynamics is identical for the 1 01 ) and 1 10) 
subspaces, so we will concentrate only on the 
levels within the 1 01 ) subspace. In the basis 
{|01), |0X_), \X+1), \X + X_)} the Hamiltonian is 



H = 





n/2 
m/2 




n/2 
s 

Vf 

m/2 



m/2 
v F 
5 

tt/2 





m/2 
n/2 

26 + V XX J 



(21) 



where Vp = 2eMih : hh- We follow the evolution of the 
four levels within this subspace over the course of the 
excitation pulse. We are interested in the populations of 
the 1 01) and |X + X_) at the end of the excitation pulse, 
and so in Fig. [4]we plot the evolution of these populations 
over time for a range of realistic mixing factors. 

We see that, as a result of the strong biexciton shift 
Vxx, the population build-up of the |X_X + ) is sup- 
pressed for a range of realistic mixing factors. For mixing 
factors upto e = 0.05 the computational states 1 01 ) and 
1 10) are also effectively depopulated - and so we conclude 
that hole mixing not a serious problem for the parity pro- 
jection. 



B. QD Spatial Separation 

The two quantum dots that form our CQD structure 
are naturally spatially separated. As has been mentioned 
in Section \\H\ any source of distinguishability will affect 
the coherence between the odd-parity states upon detec- 
tion of a photon. To analyse this effect we derive a new 
CME describing the detection step, which includes the 
effects of spatial separation. 

This is done in three stages: firstly we derive a Marko- 
vian master equation from a microscopic Hamiltonian de- 
scribing the full dynamics of the system and bath. We 
then define a jump super-operator which describes the 
evolution of the system upon a detection event. Finally 
we are able to identify a CME by imposing the condi- 
tion that the time averaged CME is equal to the Master 
Equation derived in the first step. 

The Markovian Master equation is derived from the 
following integro-differential equation: 



f 

Jo 



drTr ph {[fr/(t), [Hi{t - r),p(t) ® p ph }}} 

(22) 

where Hi(t) is the microscopic Hamiltonian in the inter- 
action picture, p(t) is the density matrix for the system 
and pph is the density matrix for the photon bath. If our 
interaction Hamiltonian is of the general form: 



(23) 
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FIG. 4: Populations of (a) 1 01) and (b) \X-X+) levels through the duration of the excitation pulse. The parameters are chosen 
as follows: Mhh,hh = Mih,hh = Vf and Ihh = hh- The laser is tuned to the exciton creation energy ooq so that 6 = 0. 



then we may write the master equation in the Born- 
Markov approximation as: 



POO 

/ drY,c a A°MUt)Mt-T)p(t) 



-A p (t-T)p{t)Al{t)} + h.c. (24) 



where C (x ^{s) = Tr p h[B^(s) B p(0) p p h] is the enviroment 
correlation function. 

For our case of two quantum dots coupled to a photon 
bath, the microscopic Hamiltonian is: 



H = H, 



CQD 



■if, 



ph 



Hi* 



with 



Hcqd 



t, 



Y2/(k)a k e ikr (4 



(25) 



(26) 



X )+H.c. 



cxo,ox represents the annihilation operator for an ex- 
citon on dot A, B respectively and is the annihilation 



operator for a quantum of the electric field, k is the 
wavevector for the electric field, /(k) = (/Iflk^k with jl 
the dipole moment vector for each qubit, <Tk the polari- 
sation vector for the electric field and ek is the energy of 
a mode k of the electric field. Finally, Ar is the center- 
to-center separation of the two quantum dots. 

We first transform to the interaction picture defined 
by Hq = Hcqd + H ph : 



Hj = ^/(k)e ik - r e i ^-°) t a k (4 + e* k - Ar cJ x ) + h.c. 

(27) 

We now proceed to calculate the master equation using 



expressions Eq. [23] and Eq. [24] We must first identify the 
system and bath operators, A{(t) and Bi(t) respectively. 
We choose: 



4(t) = /(k)e ik - r (4 + e ikAr ct x )e-o* = 

Bl(t) = a k e-^ kt . (28) 

The master equation then becomes: 



p(t) = -2ttJ2 |(1 + N(u k ))[PlP k p(t) - 2P k p(t)Pl + p(t)PlP k ]5(u;o - a*) 



-N(u; k )[P k Plp(t) - 2Plp(t)P k + p(t)P k Pl\5(LO - wk). 



(29) 



As < Wfe we may assume that N^^) = 0, and so the master equation back in the Schrodinger picture is: 

p = -i[H,p}+^(P k pPt~{PlP k p+pPtP k })5(w -u; k ). 

(30) 
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Our detection model is simply a single photon detection 
so we now can define our jump super-operator as: 



(31) 



where the fc'-vector runs over the solid angle covered by 
the detector. It is assumed for simplicity that the detec- 
tor covers the full solid-angle and also assuming an over- 
all detector inefficiency 77, the stochastic master equation 
becomes: 



k 

~{PtP k ,p}}5(LU -LO k ). 

Performing the sum over all modes, we obtain 



■i[H,p} + (l- V )Jp-Ap 



(32) 



(33) 



where 



Jp = T^cxopc^Q + 3f(k Ar)(c X opcl x + c ox pc xo ) 

+coxp4x!» ( 34 ) 
Ap = T 1 [c xo c X oP + pc xo c x + c\ x cqxP + pc\ x c ^x] 

(35) 



and 



/(«) 



2acos(a) + (a 2 — 2) sin(a) 



or 



(36) 



The function 3f(koAr) characterizes the decohering 
effect of distinguishable photons. It takes a value of 
unity for perfectly indistuishable dots, and then Eq. 33 
reduces to Eq. [7] Stacked self-assembled QDs have sep- 
arations of the order 5 nnP, while the typical exciton 
creation energy is c^o = 2 eV. This gives a value of 
3f(koAr) = 0.99925, and we can conclude that spatial 
separation has a negligible effect on the successful oper- 
ation of the parity- measurement. 



C. Detuning of QD excitonic energy levels 

We have up to this point neglected any inhomogeneity 
in the underlying structure of the two dots in our coupled 
system. However, in practice we must expect a certain 
degree of inhomogeneity: for example, due to the growth 
technique self-assembled vertically coupled quantum dots 
will tend to have different sizes. This will result in differ- 
ences in the confining potentials for the two dots, which in 
turn will impact on the exciton creation energies, and the 
overlap integrals which determine the coupling with the 
laser field. For our parity measurement, the most impor- 
tant effect will come from the detuning or non-resonance 
of the exciton creation energies for the two dots. This 
will affect both our ability to perform the excitation step 
and also the ability to retain coherence when we measure 
a photon. These effects will be analysed in this section. 



1. Excitation Pulse 

To begin, let us concentrate on the excitation step. 
As with the hole-mixing, we have two primary concerns: 
firstly we must maintain the population of the even- 
subspace within the computational ground states, and 
secondly we must ensure that the populations within the 
odd-subspace are completely transferred to the excitonic 
levels. 

Let us initially focus on the first issue. The 
Hamiltonian for the 11 subspace, written in the basis 
{|11),|1X),|X1),|XX)}, is now: 



H 



( n/2 n/2 

n/2 s B v F 

n/2 v F s A 

V n/2 n/2 s A 




(37) 



where 5a,b = va,b - ml- 

As before, we begin by transforming the Hamiltonian 
into a basis of n = eigenstates. The single exciton 
subspace is transformed using the following: 



\1X) 
\X1) 



COS ( 



— suit 
cos# 



IV--) 
!</<+> 



(38) 



where 6 = arctan(2Vp/(5 i 4 — 5b)) is the mixing an- 
gle for the two states {|1X),|X1)}, and {|^_), |^ + )} 
are the new eigenstates. Reintroducing the laser cou- 
pling, we find that the Hamiltonian in the new basis 
{|ll),|^_),|^ + ),|XX)}is: 



H = 



( n_/2 n + /2 

n_/2 s B 
n + /2 s' A 

\ n_/2 n + /2 s A 



\ 

n_/2 

n + /2 

- S B + Vxx ) 



(39) 



where Q± = ^(cos^ ± sin#) and 



)a cos 
5a sin 2 



Ob sin 

5b cos 2 



V F sin20 

V F sin 20 (40) 



To suppress any transitions to the excitonic levels, we 
require that 



n±/2 < 5 A , 5 B , 5 A + S B + V XX . 



(41) 



For both limits of small detuning (5 = 5a — 5b < 2V F ) 
and large detuning (5 = 5a — 5b > 2Vp), this condition 
is satisfied. 

Let us now turn our attention to the dynamics within 
the odd-subspace during the excitation phase. For a 
successful outcome, we must effect the same population 
transfer into the exciton levels for both the 1 01 ) and 1 10) 
initial states. Our Hamiltonian for this subspace in the 
rotating frame is: 



H = 



( n/2 

n/2 5 A 

ft/2 

V on/2 5 B 



(42) 
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2. Detection 
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FIG. 5: Plot of the population of the single excitonic levels 
and the phase between them at the end of the excitation step 
as a function of the detuning of the two QDs 



We next discuss the effect of detuning on the detec- 
tion process. To analyse this we will derive a master 
equation from first principles describing the relaxation 
process. From here we will proceed to identify the rel- 
evant jump operator associated with the detection of a 
photon, and the resulting CME governing the dynamics 
of the system until the measurement occurs. 

We start with a microscopic Hamiltonian describing 
the coupled quantum dots, the photon bath, and also 
the coupling between these two systems: 



H = H, 



CQD 



ph 



■H, 



int 



(43) 



where 



Hcqd = uac^qCxo+ubcIxCox 

H ph = ^ UJ k a L a * 
k 

Him = ^2g^al(c ox +c X o) +H.c. (44) 



To achieve equal population transfer to the |0X) and 
|X0) the laser is tuned to the midpoint between the two 
excitonic levels: 5 a = —5b = 5 = (uoa — ^#)/2. The 
maximum population of the excitonic levels is given by 
ft 2 /(ft 2 +c5 2 ) and is found at a time ofr,, = + 5 2 . 

For Q ^> c5 we can achieve the required tt pulse for each 
odd-parity state. However, we must also ensure that the 
coherence between the 1 01) and 1 10) states is preserved 
as we transfer the population to the excited states. Since 
the levels are not resonant, a net phase difference can be 
accumulated between them by the end of this excitation 
pulse. The variation of this phase is shown, along with 
the population of the excitonic states, in Fig. [5] as a 
function of detuning. The figure seems to show that we 
require very similar dots to ensure a successful excitation 
pulse. However, if the phase accumulated is known, it is 
possible to correct for at the end of the measurement 
using single qubit operations. 



where g(k) is the photon-exciton coupling constant. Note 



that in contrast to Eq. [27] , here we are neglecting spatial 
separation but accounting for frequency discrepancy. 

We begin by transforming into the interaction picture, 
defined by H = Hcqd +H p h. The resulting interaction 
Hamiltonian is: 



Hi = 5^0k4( 



t ( e i(^-u A )t Cxo + e i(^-u B )t Cox ^ + h c 



(45) 



We proceed as before (Eqs. |23|and 24) by defining the 
system and environment eigenoperators Ai(t) and Bi(t), 
respectively: 



4(t) = //(kH<- • , 



xo 



), 



K(t) = a k e-^ k£ . (46) 



Assuming that we are operating at zero temperature, we 
obtain the following master equation: 



p(t) = r A c xo p(t)c xo + (T A + r B )/2{c X op(t)4x^ UB ~ u,A)t + c xp(t)4 e-^-^) 

+r B coxp(£)4x - ^{c xo c X o + 4x c ox, pit)} (47) 

I 



where the decay rates of the two dots are given by 

T Ay B = 27r^ k |#k| 2 <K^k - va,b)- If {u B ~ u A ) > 
(r^ + r#)/2 the fast oscillating terms may be neglected 
as their contribution would average to zero on the time 



scale of the relaxation. In this case, all coherence is 
lost during the monitored relaxation. Therefore, for a 
non-destructive measurement, we must work in a regime 
where (uob— u>a) <^ (r A+r /2. This condition is equiv- 
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alent to requiring that there is a large overlap of the two 
spectral lines from the two quantum dots, so that the two 
dots are to a high degree spectrally indistinguishable. 

We may make a further assumption to simplify the 
analysis: the decay rates Ta,b are proportional to uo\ B , 
and so if the detuning is small we may assume that the 
decay rates are identical for the two dots, ie Ta = T# = 
T. The resulting master equation back in the Schrodinger 
picture is then: 



p(t) = -i[H C QD,p] + r(c x + cxo) P (clx + c 

r 



XOJ 



(48) 



To derive the conditional master equation, we identify 
the jump operator associated with a photon detection 
event. As we are assuming that our detector covers the 
full 47r solid angle, our jump operator is simply J[p(t)] = 

T(cox + cxo) P (clx + c xo)- Assuming some loss in the 
detection process, again parameterised by the variable 77, 
the unnormalised CME is 



p = -i[H CQD ,p] + (1 - rf)Jp - Ap 



(49) 



where Ap = |{ c xo c ^o + cj x c x, p}- 

We see that when a photon is detected at a time t^, 
a phase of e*^ - ^^ is introduced between the 1 01) 
and 1 10) states. This phase is the relative phase accumu- 
lated during the time that the system spends in the exci- 
tonic levels. This is a general problem exhibited by many 
optical-matter measurement based schemes. Although, 
in theory we could have information about the relative 
phase introduced between the two levels, in practice this 
information can beyond our reach, and effectively intro- 
duces a random phase. Such a random phase will de- 
stroy any coherence between the two states in the odd- 
subspace, thereby destroying the non-destructive parity 
measurement. 

In order to retain coherence during the measurement, 
we therefore require a photon detector with good enough 
time resolution that we are able to successfully access 
this phase information. If this is possible, then we may 
correct for the phase accumulated using single qubit ro- 
tations, for example with a single Z-rotation on the first 
QD. 

In the worst-case scenario of distinguishable dots, we 
have a probabilistic scheme for generating entanglement, 
which is successful only when we project into the even 
subspace. However, this still provides a powerful re- 
source which may be used for scalable QIP. Furthermore, 
in this limit of distinguishable dots, the scheme is compa- 
rable with many of the existing, inherently probabilistic, 
schemes for optical-matter measurement-based QIP^. 



VI. EXPERIMENTS 

We will conclude by presenting a discussion of exper- 
imental procedures for testing the fidelity of the parity 
measurement. 



Quantum process tomography provides a general pro- 
cedure for characterising the dynamics of a quantum sys- 
tem provided that we can measure each qubit indepen- 
dantly in the X, Y and Z bases. However, in our system, 
we must restrict ourselves to only measuring both qubits 
in the same basis at the same time, since the two dots 
are not each individually addressable, and we may only 
perform global single qubit rotations. 

We first introduce entanglement witnesses for the four 
possible entangled Bell states. It has recently been 
shown that entanglement witnesses can be constructed 
for highly entangled states using the stabilisers that de- 
fine these states^. An observable Sk is a stabiliser for 
the state 1-0) if: 



W> = \i>)- 



The stabilisers for the four Bell states are then: 



(50) 



Bell States 



Stabilisers Sk 



KM = ^|oi) 
1 |oi) 



!</>+> 
10-) 



± 



1 100) 
1 100) 



|10) 
|10) 

111) 
111) 



— Z1Z2, Xlx 2 
-Z1Z2, —XIX2 
Z1Z2, XIX2 
Z\Z2-, —X1X2 



By making measurements on XX and ZZ we are able 
to distinguish between the four Bell states. The ZZ mea- 
surement is related to the parity projection operators by 
Pe,o — I±z ^ z<2 . Meanwhile the XX measurement is 
achieved by first performing a Hadamard gate on each 
qubit, and then performing the same spin-parity mea- 
surement. This XX measurement effectively measures 
the phase between the states within each subspace. We 
therefore are using repeated applications of the parity 
measurement to gain information about sits own opera- 
tion. 

As we have seen from the previous analysis (section[V|), 
the primary error source in this system is decoherence be- 
tween the states in the odd-subspace due to distinguisha- 
bility of the dots. We can therefore assume that there is 
some loss of coherence after the parity measurement, as 
follows: 



Popph 

PepPI 



PoipP^ + aPoipPio + aPwpPoi + PwpPu 
(Poo + Pu)p(P^ + p!i)- 



10' 



(51) 



where Pi = a denotes the degree of coherence; it 

takes a value of unity for indistinguishable dots, and will 
be somewhat less than that for distinguishable dots. 

Starting with an initial state that is an equal superpo- 
sition of the four computational states, the state of the 
system conditioned on observing a photon is: 



Pi 



/ 

1 a 

a 1 

\ 



(52) 
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Then after a the global Hadamard rotation, the state 

is: 

/ 1 + a -1 - a \ 
1 1 — a a — 1 /roN 

" I= 4 a-1 l-a ' ^ 53) 
\ -1 - a l + a J 

Finally, the probability of projecting into the odd- 
subspace after the second parity measurement is 
We therefore see that any loss of coherence during the 
parity measurements will manifest itself in the final prob- 
ability, thus giving us a clear method of quantify the ef- 
fect of distinguishability on the non-destructive nature of 
the parity measurement. 



coupled quantum dots. We have estimated the fidelity 
of the parity measurement scheme presented here in the 
presence of realistic sources of errors. We find that the 
measurement is robust in the presence of inefficient de- 
tectors, ineffective spin-selective excitation and spatial 
separation of the dots. For spectrally separated dots, it 
is found that the performance of the measurement is de- 
pendant on the degree of overlap of the spectral lines from 
the two dots. Total spectral distinguishablility results in 
a probabilistic measurement which is still sufficient for 
growing large scale entangled states. Finally, we have 
proposed an experimental method that is able to verify 
the success of the parity measurement and quantify the 
degree to which the measurement can be performed in a 
non-destructive manner. 



VII. CONCLUSIONS 



In conclusion, we have presented a novel scheme for 
implementing a spin-parity measurement on a pair of 
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